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Experimental  Area 

CN2-optics  vs  Cj^2-bulk  using  bulk  measured  sea 
surface  temperature.  Solid  line  is  perfect  agreement, 
dashed  lines  are  factor  of  two  disagreement. 

Aircraft  measured  IR  sea  surface  temperature  and 
Cip2  vs  distance  from  shore.  Optical  path  location 
is  shown  by  a  vertical  dark  line. 

CN2-optics  vs  Cjf2-turbulence,  data  not  correct¬ 
ed  for  salt  Incrustation.  Solid  line  is  perfect 
agreement,  dashed  lines  are  a  factor  of  two 
disagreement . 

CN2-optics  vs  C^2-turbulence ,  data  corrected 
for  salt  encrustation.  Solid  line  is  perfect 
agreement,  dashed  lines  are  factor  of  two  disagreement. 

CN2-optics  vs  Cj»j2-bulk  using  IR  measured  sea 
surface  temperature.  Solid  line  is  perfect  agreement, 
dashed  lines  are  factor  of  two  disagreement. 

Difference  In  measured  sea  surface  temperatures 
(IR-bulk)  vs  time.  Shaded  areas  indicate  cloud  cover. 
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I.  INTRODUCTION 

Surface  layer  turbulence  models  have  been  improved  to  the 
point  where  quite  good  estimates  of  turbulence  intensities  can  be 
made  from  readily  measured  meteorological  parameters  (wind, 
temperature,  and  humidity).  However,  applications  of  present 
formulations  to  estimate  turbulent  intensities  of  the  optical 
index  of  refractions  within  the  surface  layer  have  been  limited  by 
acknowledged  deficiencies.  These  are  the 

1.  bulk  scaling  of  the  contributions  of  the  turbulent 
variance  of  water  vapor  and  the  turbulent  covariance  of 
water  vapor  and  temperature. 

2.  observational  verification  of  the  role  of  turbulence  on 
overwater  optical  degradation  based  on  both  optical 
measurements  and  model  estimates  of  Cpj2,  the  refractive 
index  structure  function  parameter. 

The  purpose  of  this  report  is  to  describe  combined  overwater 
measurements  that  allow  a  comparison  of  optical  measured  values  of 

with  values  calculated  from  both  turbulence  and  bulk 
measurements.  The  turbulence  estimates  include  contributions  from 
both  variance  (C-j2  and  Cq2)  and  covariance  (Cijq)  components. 

The  bulk  model  is  that  formulated  by  the  Naval  Postgraduate 
School  (NPS)  and  utilizes  stability  corrected  scaling  parameters. 
As  will  be  shown  here  the  bulk  formulation  yields  results  which 
are  considerably  better  than  direct  turbulence  measurements  and 
are  in  good  agreement  with  optical  measurements. 
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In  order*  to  meet  the  above  needs  the  Naval  Environmental 
Prediction  Research  Facility  and  NPS  planned  a  series  of  coinci¬ 
dent  optical  and  meteorological  measurements  to  be  made  on 
Monterey  Bay.  This  work  was  undertaken  during  the  Marine  Aerosol 
Generation  and  Transport  (MAGAT)  experiment.  The  work  was  per¬ 
formed  from  28  April  to  9  May  1980  by  the  Environmental  Physics 
Group  of  NPS  in  cooperation  with  the  NPS  Optical  Propagation  Group 
and  Airborne  Research  Associates  (ARA).  The  purpose  of  the  effort 
was  to  verify  overwater  optical  propagation  models  that  have  been 
developed  to  predict  extinction  and  scintillation.  The  purpose  of 
this  report  is  to  evaluate  the  bulk  aerodynamic  method  for  obtain¬ 
ing  the  index  of  refraction  structure  function  from  mean  meteoro¬ 
logical  parameters.  Evaluation  of  the  NPS  boundary  layer  aerosol 
model  will  be  the  subject  of  another  report. 

During  MAGAT,  the  full  range  of  meteorological  measurements 
were  made  on  the  RV/ACANIA  and  on  the  ARA  aircraft.  This  included 
both  mean  and  fluctuating  parameters.  All  model  evaluations  in 
this  report  were  made  using  the  shipboard  data.  Optical  measure¬ 
ments  were  made  on  the  13  km  overwater  range,  at  regular  periods 
around  the  clock  for  eight  days  in  order  to  experience  as  wide  a 
range  of  conditions  as  possible.  The  ship  was  stationed  on  the 
optical  path  frequently  for  direct  comparison  with  the  optical 
measurements . 


6 


II.  EXPERIMENTAL  ARRANGEMENT 

The  RV/ACANIA  was  equipped  with  a  multi-level  measurement 
system  to  measure  both  mean  and  fluctuating  meteorological 
parameters.  The  heights  at  which  sensors  were  placed  above  mean 
sea  level  and  the  quantities  measured  are  shown  in  Table  1: 


TABLE  1.  Shipboard  Measurement  Configuration 


Height 

0 


4.2  m 


7.0  m 


19*6  m 


Parameters  Measured 

Sea  surface  IR  temperature  (Tjp) 

Bulk  water  temperature  (Tg) 

Mean  temperature  (T) 

Mean  wind  speed  (U) 

Wind  speed  fluctuation  ( U ' ) 

Mean  temperature 
Mean  wind  speed 
Mean  dew  point  (Tp) 

Temperature  fluctuation  ( T ’ ) 

Wind  speed  fluctuation 
Humidity  fluctuation  ( H ' ) 

Mean  temperature 

Mean  wind  speed 

Mean  wind  direction  (WD) 

Mean  humidity 

Wind  speed  fluctuation 
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In  addition  to  these  sensors,  the  visibility  (V)  and  inversion 
height  (Z^)  were  determined  and  aerosol  spectra  were  measured  at 
a  height  of  8.5  ra. 

The  sensors  used  were: 

T jpj  Barnes  PRT-5 

Tg,  T  Rosemount  platinum  thermometers  mounted 

in  Gill  aspirators 
U,  WD  MR I  1022  system 

Tj)  General  Eastern  cooled  mirror 

T'  2 . 5 y  platinum  microthermal  sensors  and 

Sylvania  140  bridge 

U'  60 p  platinum  on  quartz  substrate  and  TSI 

1054  bridge 

H'  ERC  Lyman-Alpha 

Zj  Aerovironment  300  Sounder 

V  MRI  1580  Fog  Visicmeter 

Two  T'  sensors,  placed  a  distance  of  30  cm  apart,  were  the 
primary  sensors  for  determining  0^2.  a  single  T'  sensor  was 
placed  Immediately  adjacent  to  the  active  volume  of  the  Lyman- 
Alpha  in  order  to  measure  the  temperature-humidity  cospectrum. 

The  temperature  measurement  circuitry  in  the  dew  point  sensor  was 
not  used  due  to  Instability  problems.  This  system  has  a  platinum 
sensor  mounted  in  a  three  wire  configuration.  The  leads  were 
changed  to  4  wire  to  improve  the  accuracy  and  the  same  system  was 
then  used  to  measure  all  of  the  temperatures  Tg,  T  and  Tp. 


The  IR  thermometer*  was  mounted  on  a  railing  on  the  ship  approxi¬ 
mately  4  ra  above  the  water.  The  sensor  was  angled  at  approximately 
45°  to  insure  that  the  ship  wake  was  not  included  in  the  field  of 
view.  The  platinum  thermometer  was  inserted  in  a  brass  plug  in  the 
end  of  a  long  1"  diameter  tygon  tube.  The  arrangement  was  slightly 
less  buoyant  than  desired  and  floated  so  that  it  averaged  the  water 
temperature  for  about  the  first  12"  below  the  surface.  The  depth 
depended  on  the  ship  speed. 

All  data  were  recorded  with  a  Hewlett  Packard  3052  data  acqui¬ 
sition  system  controlled  by  a  Hewlett  Packard  9825S  computer. 

Almost  all  data  were  obtained  as  one  half  hour  averages.  The  same 
voltmeter  was  used  to  measure  voltages  and  the  4-wire  resistances. 

Neither  descriptions  of  the  optical  measurements  system^) 
nor  presentation  of  the  resulting  data^2)  will  be  given  here 
as  they  are  included  in  another  report.  Figure  1  shows  the  location 
of  the  measurement  area  in  Monterey  Bay.  The  optical  path  is 
located  so  that  it  is  approximately  perpendicular  to  the  prevailing 
northwest  flow  in  such  a  location  that  land  influence  is  minimal. 
Also,  scintillation  measurements  weight  the  center  of  the  path 
further  reducing  any  land  influence.  For  most  of  the  measurements 
reported  here  the  ship  was  anchored  at  the  position  shown  in  the 
figure;  a  few  measurements  were  made  while  the  ship  was  in  motion 
within  the  square  area  shown.  When  anchored,  the  ship  automatically 
faced  into  the  wind  and  while  underway,  data  was  only  taken  when  the 
ship  was  headed  into  the  wind. 


The  ship  system  measures  the  properties  of  the  air  that  passes 
its  location  during  a  one  half  hour  period.  The  optical  path 
averages  the  properties  of  the  air  over  the  13  km  path  length.  We 
are  thus  comparing  a  time  average  and  a  space  average.  Or,  if 
there  is  horizontal  homogeneity  along  the  mean  wind  direction,  the 
time  average  taken  at  the  ship  is  equivalent  to  a  space  average 
along  the  mean  wind.  In  any  event,  the  two  averages  obtained  by 
the  optical  technique  and  by  the  meteorlogical  measurements  are 
not  exactly  equivalent.  Exactly  how  this  affects  comparisons  of 
optical  and  meteorlogical  results  cannot  be  determined  without 
detailed  knowledge  of  the  local  airflow  and  temperature  patterns. 
This  topic  is  now  under  investigation. 

Table  1  lists  the  meteorological  data  that  was  acquired 
coincident  with  optical  scintillation  measurements.  Listed  are 
the  true  wind  speed,  sea  temperature  from  both  bulk  and  IR 
thermometers,  air  temperature,  relative  humidity,  and  the 
calculated  stability  parameter  Z/L. 


11 


oY»  T  0S4T-8Q/S0 


III.  BULK  MODEL 


The  bulk  model  uses  the  differences  In  the  values  of  mean 
parameters  between  the  surface  and  a  reference  height  to  estimate 
small  scale  properties  of  the  atmosphere.  This  MPS  bulk  model  was 
first  developed  In  1977  to  calculate  the  Index  of  refraction  struc¬ 
ture  function,  C^2,  and  was  first  applied  to  evaluate  results 
from  the  CEWCOM-78  experiment^).  Since  that  time,  NPS  has 
verified  the  overwater  scaling  and  stability  correction  functions 
for  wind,  temperature,  and  humidity.  The  model  has  proven  to  work 
quite  well  for  predicting  small  scale  fluctions  in  wind,  temperature 
and  water  vapor  as  was  shown  by  comparing  direct  fluctuation 
measurements  with  the  bulk  calculated  values . ( ^ » 5, 6) 

The  CEWCOM-78  report  contains  a  sketch  of  the  model  and  the  full 
development  Is  presented  below. 

The  optical  refractive-index  structure  function  parameter  is 
related  to  the  temperature  structure  function  parameter,  C^2, 
and  the  humidity  structure  function  parameter,  Cq2,  by  (7) 

CN2  =  (79  x  10-6  P/T2)2  (CT2  +  .113  CTQ  +  3.1  x  10~3  CQ2)  (1) 

Where  P  is  the  pressure  in  mb  and  T  is  the  absolute  temperature. 

CTq  is  the  temperature-humidity  cospectral  structure  function. 

A.  Monin-Obukhov  Scaling 

Cx2  can  be  related  to  the  measured  meteorological  quantities 
through  Monin-Obukhov  surface  layer  similarity  parameters^, 9) 

CT2  =  T#2  Z-2/3  f(£), 

Cq2  -  Q#2  Z-2/3  a  f(g). 


(2a) 

(2b) 
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where  T#  is  the  potential  temperature  scaling  parameter,  Q*  is 
the  water  vapor  density  scaling  parametr  (gn/m3),  z  is  the 
height  above  the  surface,  £*  Z/L  is  the  similarity  ( dimensionless ) 
height  parameter  and  f( O  is  the  empirical  function  found  by 
Wyngaard,  et  al.  The  quantity  A  is  a  constant  approximately  equal 
to  0.8. (5)  The  cospectral  function  is  given  by 

C?q  *  r*TQ  T#Q#  Z”2/3  a^2  f(£).  (2c) 

where  r^Q  is  the  temperature-humidity  correlation  parameter 
equal  to  0.8  under  unstable  conditions.  The  value  of  r-pg  for 
stable  conditions  is  not  well  known  in  the  surface  layer.  Note 
that  Q#  in  (gm/m3)  and  q*  in  (gm/kg)  are  related  by  Q*  *  1.3  q* 
at  the  surface,  q#  is  the  water  vapor  mixing  ratio  scaling 
parameter.  The  Monin-Obukhov  length  scale,  L,  is  defined  by 

L  a  (TAg)  U*2[T#  +  6.1  x  10 _1*  Tq*]-1,  (3) 

where  k  is  von  Karman's  constant  (0.35),  g  is  the  acceleration  of 
gravity,  and  U#  is  the  friction  velocity. 

The  problem  of  predicting  Cj^2  is  now  reduced  to  finding 
values  for  q#,  T*  and  £  (or  L).  The  bulk  method  is  based  on 
obtaining  values  of  temperature,  relative  humidity,  and  wind  speed 
at  the  sea  surface  and  at  some  reference  height,  Z.  The  differ¬ 
ence  between  the  surface  value  and  the  value  at  height  Z  can  be 
related  to  the  scaling  parameter  through  the  profile  equations . ( 10 ) 
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U#  *  kU  l in  Z/Z o  C)  ]_1, 


T#  -  (T  -  Ts)  aT  k  lln  Z/Z0T  -C2  ( C)  ]~1 


Q#  =  (q  -  qs)  ci ip  k  [£n  Z/Zqt  ~C  2^)3”^  (4c). 

where  >p  in  the  ratio  of  heat  transfer  to  momentum  transfer  at  £  *  0, 
and  5  is  the  value  at  height  Businger,  et  al.  found 

cup  3  1.35,  others  have  found  different  values.  The  quantities  Z0 
and  Zgip  are  the  roughness  lengths  for  velocity  and  temperature 
profiles.  Mote  that  these  equations  can  be  written  in  the  standard 
drag  coefficient  form 


U*  3  cqI/2  u> 


'*  =  ct1/2  (T  -  Ts), 


q *  ■  c-p1/2  (q  -  qs). 


In  Equs .  4  and  5,  we  have  assumed  that  the  water  vapor  dependences 
(q)  can  be  treated  with  the  same  coefficients  as  the  temperature 

(  Zgip,  Cip)  . 

The  stability  dependence  of  the  drag  coefficients  can  be 
obtained  from  Equs.  4  and  5 


cDl/2  =  (k/£n  Z/Z0)[l  -  (£n  Z/Zq)"1  ^CO]"1 
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CT1/2  =  <aT  k/£n  Z/Z0T)[1  -  <£n  Z/Z^)"1  i^U)]-1 


(6b) 


We  can  define  the  neutral  stability  drag  coefficients  in  terms  of 
the  roughness  lengths  as 

(7a) 

(7b) 

Note  that  the  given  drag  coefficient  at  height  Z,  one  can 
calculate  the  roughness  length 


cdn1/2  =  k(£. n  Z/Z 0)"1, 
Ctn^^2  *  aT  ^  Z/Zqt)"1 


Z0  -  Z  e  xp[  -  k/c^M1/2] 
Zqiji  =  Z  exp[-oiip  Ic/ctn^^2 


We  are  now  able  to  calculate  the  atmospheric  stability  at 
height  Z,  5  =  Z/L,  using  Equs.  3,  4  and  7, 

5  =  C0  [1  -  Un  Z/Zq)'1  ^  (C)]2[l  -  Un  Z/Z^)*1  *2U)]  (9) 

where 

C0  =  (kgS/T)  (ctn1/2/<=dn)  (AT  +  6.1  x  10~4  TAq)u"2  (10) 
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B.  Empirical  Constants  and  Quantities 

We  have  been  using  a  value  of  von  Karman's  constant  k  ■  0.35 
based  on  the  original  Businger,  et  al .  work.  Recently, 

Garratt^^)  has  published  a  survey  which  implies  k  =  0.41. 
Businger,  et  al.dD  found  a<p  =  1.35,  however,  if  one 
uses  k  *  0.41  then  a  value  of  aT  =  1.15  would  maintain  a 
constant  onp  k  . 

A  typical  value  of  cDN  at  Z  3  10  m  is  1.3  x  10“3  which 
yields  Zg  3  6  x  10-i<m.  Kondod-3)  and  Garratt^^) 

both  give  equations  for  wind  speed  dependence  of  the  Z  =  10  m  drag 
coefficient.  Kondo's  formulae  are  used  in  our  model  formulation 
and  are  given  in  Table  1. 


Table  2.  qjN  versus  wind  speed  at 
UCrns-1) 

.3  -  2.2 
2.2  -  5.0 
5.0  -  8.0 
8.0  -  25.0 


10  m  from  Kondo^1^). 

CDN  x  10 ^ 

1.08  x  U--x5 

.77  +  .086  x  U 
.87  +  .067  x  U 
1.2  +  .025  x  U 


The  temperature  drag  coefficient  has  been  measured  by  several 
groups  (see  Davidson,  et  al.(^),  for  a  summary),  but  we  feel 
a  best  estimate  is  cp^  *  1.3  x  10“3  at  Z  =  10  m.  Assuming 
am  *  1.35,  we  obtain  Zqt  =  2.0  x  10-5  m.  For  our  bulk 
model,  we  assume  that  Zqt  is  independent  of  wind  speed  and 
that  the  wind  speed  dependence  of  Zq  can  be  obtained  from 
Rondo’s  CpN  using  Equ.  8a  with  Z  *  10  m. 
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C.  Procedure 


1.  Input  data  are  sea  surface  temperature  (Tg),  air 
temperature  (T),  relative  humidity  or  dew  point  (H  or  Tp)  and 
wind  speed  (U).  The  last  three  are  measured  at  a  reference  height 
Z.  From  T  and  H  (or  Tp)  calculate  q.  From  Tg  calculate  qg 
assuming  that  H  =  100%  at  the  surface. 

2.  From  U,  calculate  (Kondo)  for  Z  =  10  m 

From  cqn>  Z  =  10,  calculate  Zq  (Equ.  8a) 

Let  Zq<p  =  2.0  x  10-5 

Let  ctn  =  1-3  x  10“3  if  Z  =  10  m.  If  Z  f  10  m, 
use  Equs .  7a  and  7b  to  calculate  the  drag 
coefficients . 

3.  From  AT  =  T  -  Tg  (potential  temperature) 

Aq  =  q  -  qg 

AU  =  U,  calculate  (Equ*  1°) 

4.  Solve  Equ.  9  Iteratively  to  obtain  Z  from  £o*  Note 

that  L  =  Z/5 • 

5.  From  T»,  Q*  =  1.3  q#>  and  Z/L  calculate  0^- , 

Cq^  and  C<tq  at  any  height  using  Equ.  2. 

Calculate  from  Equ.  1. 


D.  Stability  Correction  Functions 


Velocity  Profile: 

4^(5)  -  2  in  [(1  +  x)/2  ]  +  in  [(1  +  x2)/2] 

-2  tan  -1  (  x)  +  tt/2  C<0 

x  =  (1  -  15  5  )1/4 

4»l(5)  -  -JU7  5  5>0 


Temperature  Profile: 

4>2(C)  =  2^n  [{1  +  x)/2]  C<0 

x  «  (1  -  9£)1/2 

♦  2(5)  =  -6.55  5>0 
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IV.  COMPARISON  WITH  OPTICAL  MEASUREMENTS 

Table  3  presents  the  following  computed  quantities:  (1) 
scaling  parameters  for  wind,  temperature  and  water  vapor,  (2) 

Cpj2  from  turbulence  measurements,  (3)  C  from  the  bulk 
model,  (4)  the  optically  measured  Cj^  and,  for  reference,  the 
stability  parameter. 

Comparisons  of  calculated  from  the  bulk  model  with 

those  measured  optically  are  shown  in  Figure  2.  The  solid  points 
are  for  cases  where  the  surface  layer  is  in  equilbrium  and  the 
open  circles  are  for  non-equilibrium,  which  will  be  explained 
below.  The  solid  line  indicates  perfect  agreement  and  the  two 
dashed  lines  are  for  a  factor  of  two  disagreement.  For  all  but 
two  of  the  eighteen  cases  where  the  surface  layer  was  in 
equilibrium  the  agreement  is  within  a  factor  of  two.  The  mean 
percent  error,  taking  the  optical  value  to  be  correct,  for  all 
equilibrium  values  is  33$.  This  is  very  good  agreement. 

We  have  found  that  at  times,  there  is  a  change  in  water 
temperature  in  Monterey  bay  in  the  vicinity  of  the  optical  beam. 
The  change  is  from  colder  to  warmer  as  the  shore  is  approached. 

As  implied  above,  the  change  in  temperature  is  not  always  present 
and  the  frequency  of  occurence  has  not  been  determined.  Figure  3 
shows  a  plot  of  sea  surface  temperature  vs  position  measured  by 
the  ARA  aircraft  using  an  IR  thermometer.  The  aircraft  flew  a 
course  perpendicular  to  the  optical  path  from  the  shore  to  25  km 
at  sea.  The  location  of  the  optical  path  is  shown  in  the  figure 
by  a  heavy  vertical  line.  The  water  temperature  is  seen  to 
gradually  increase  (but  not  monotonically )  by  a  few  tenths  of  a 
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Table  3.  Calculated  Scaling  Parameters  and  Index  of  Refraction  Structure  Function 


degree  up  to  about  6  km  from  shore  then  increases  by  one  degree  in 
4  km.  The  optical  path  is  in  the  middle  of  the  rapid  change 
region  on  the  day  the  profile  shown  was  determined. 

It  was  fairly  easy  to  determine  the  days  when  the  temperature 
discontinuity  existed  from  shipboard  measurements.  The  shipboard 
operation  required  that  it  move  in  and  out  of  the  bay  frequently, 
so  we  were  able  to  compare  bay  temperatures  to  those  further  at 
sea.  On  several  days  the  bay  temperature  at  the  optical  path  was 
approximately  one  degree  elevated.  We  assume  that  the  surface 
layer  may  not  be  in  equilibrium  at  the  ship  site  when  the 
temperature  jump  exists  and  have  indicated  data  from  such  days  by 
open  circles  in  Figure  2.  Note  that  during  such  times  the  optical 
path  may  also  be  inhomogeneous. 

The  results  clearly  demonstrate  the  air  is  not  in  equilibrium 
with  the  sea  surface  temeprature  at  the  ship's  location. 

Calculated  values  of  CN2  average  a  factor  of  4  times  the 
optically  measured  values.  This  is  completely  different  than  the 
results  for  the  equilibrium  cases.  The  explanation  is  as  follows: 
at  the  ship  the  elevated  sea  surface  temperature  is  measured  which 
results  in  a  large  calculated  C^2  due  to  the  large  air-sea 
temperature  difference.  However,  the  thermal  turbulence  in  the 
air  is  significantly  influenced  by  the  temperature  difference 
further  upwind.  How  large  the  calculation  will  be  in  error 
depends  on  the  time  the  surface  layer  has  to  adjust  to  the  new 
temperature. 

Further  supporting  evidence  for  this  effect  was  obtained  by 
the  aircraft  measurements  of  C?2 ,  (Figure  3).  From  25  km  at 


sea  to  the  region  of  the  optical  beam  Crp2  is  fairly  constant 
then  rising  about  70%  closer  to  shore.  The  thermal  turbulence 
does  not  respond  instantly  to  the  temperature  change,  as  was 
suggested  in  the  preceeding  paragraph. 

The  final  conclusion  concerning  the  bulk  model  is  that  it 
works  quite  well  for  predicting  optical  scintillation  for  an 
equilibrium  surface  layer.  In  non-equilibrium  situations,  the 
calculation  can  be  expected  to  be  in  error,  with  the  maximum 
error  depending  on  the  magnitude  of  local  mean  parameter 
discontinuities.  In  the  open  ocean,  where  surface  temperatures 
tend  to  be  horizontally  uniform,  the  non-equilibrium  situations 
are  expected  to  be  uncommon. 
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V.  THERMAL  TURBULENCE  MEASUREMENTS 

As  was  indicated  in  Section  II,  measurements  of  thermal,  wind 
speed,  and  water  vapor  turbulence  were  made  during  MAGAT.  These 
data  are  used  to  calculate  the  scaling  parameters  T*,  U*,  and 
q*  and  the  structure  functions  C^2,  Cq2,  and  Cu^.  These 
parameters  are  related  through  Equs  2  for  T  and  q  and  similarly 
for  U: 


Cy2  =  U,2  Z-2/3  gU). 


(11) 


Traditionally  one  uses  the  rate  of  velocity  turbulence 
dissipation,  £  ,  rather  than  Cy2,  and  they  are  related  by: 

C u2  *  2. 0  e  2/3.  ( 12 ) 


The  dissipation  stability  function,  $  (£),  is  introduced  to 
directly  relate  U#  and  e  as(^) 


e 


<M£) 


(13) 


Turbulence  signals  are  analyzed  in  two  ways:  (1)  spectral 
analysis  and  (2)  obtaining  the  rms  of  spatially  or  frequency 
filtered  signals.  The  spectral  method  is  based  upon  the 
assumption  of  the  "local  Isotropy"  and  the  Kolomogorov  -5/3  slope 
of  the  one-dimensional  power  spectral  density,  Fx(k) 

Fx(k)  =  0.25  Cx  k-5/3,  (14) 
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where  k  Is  the  wavenumber  and  x  refers  to  T,  U,  or  q.  Performing 
a  Fourier  spectrum  analysis  in  the  frequency  domain  (f)  and  using 
Taylor's  hypothesis  gives: 

Cx  =  4  (lf  )2/3  f5/3  PxCf).  (15) 

Using  two  sensors  spaced  a  distance  d  apart,  the  structure 
function  can  be  found  by  measuring  the  variance  of  the  difference 
in  x 


C  2  =  [x(r)  -  x(r  +  d)]2/d2//3, 

X 


(16) 


If  frequency  filtering  rather  than  the  spatial  filtering  is  used 
with  upper  and  lower  frequency  limits,  fu,  f^,  then 


f  and  k  are  again  related  through  Taylor's  hypothesis.  Using 
Equ.  15,  the  structure  function  Is  related  to  the  rms  signal  by 


(ll) 
v  u  ' 


2/3 


(x'  )* 

rms 

<v2/3  -  *7™* 


(18) 


These  analyses  only  apply  in  the  inertial  subrange  so  that  the 
probe  spacing  for  the  spatial  filtering  technique  and  the 
frequency  band  for  the  frequency  filtering  technique  must  insure 
that  only  this  range  Is  included. 
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Measurements  of  C<r2  by  microthermal  sensors  are  very  diffi¬ 
cult  because  of  the  problem  of  salt  loading.  This  is  due  to  the 
wires  becoming  sensitive  to  humidity  fluctuations  when  they  are  salt 
encrusted  (14)(15).  Humidity  fluctuations  will  falsely  elevate 
Crp2  and,  hence,  the  calculated  C^2.  Values  of  C^2  calculated 
from  the  turbulence  results  are  listed  in  Table  3  and  plotted  vs 
C^-optical  in  Figures  4  and  5* 

Figure  4  shows  results  for  which  no  correction  for  the  salt 
loading  effect  has  been  made  and  the  comparison  is  very  poor.  In 
Figure  5,  we  show  results  where  a  correction  has  been  made.  We 
correct  the  data  by  using  signals  that  occurred  immediately  after 
washing  the  wires.  The  comparison  improves  somewhat  but  is  still 
poor. 

If  the  thermal  turbulence  measurements  could  be  made  correct¬ 
ly,  this  method  should  be  superior  to  the  bulk  model  calculation 
since  the  small  scale  turbulence  which  is  responsible  for  optical 
scintillation  is  being  measured  directly.  However,  the  measure¬ 
ments  are  very  difficult  and  subject  to  error.  We  do  not  believe 
that  the  technique  can  be  made  to  perform  as  well  as  the  bulk 
method . 


30 


VI .  BULK  MODEL  WITH  IR  SEA  SURFACE  TEMPERATURE 

Figure  6  shows  a  comparison  of  optically  measured  and 
values  calculated  from  the  bulk  model  using  the  IR  sea  surface 
temperature.  These  results  are  not  presented  in  tabular  form. 
Obviously,  there  is  very  little  agreement  between  measured  and 
calculated  values.  (Three  points  are  off  the  graph  and  not 
plotted. )  This  means  that  the  bulk  water  temperature  is  a  better 
measure  of  the  surface  temperature  boundary  condition  for  calcula¬ 
ting  surface  layer  fluxes. 

This  is  a  surprising  result  since  the  skin  temperature,  which 
is  the  parameter  measured  by  an  IR  thermometer,  should  be  the 
desired  boundary  condition.  Note  that  in  the  bulk  model  the  bulk 
calculation,  including  the  stability,  uses  the  IR  temperature  so 
that  it  is  self  consistent. 

We  have  compared  the  IR  and  bulk  temperature  directly  to  see 
if  there  is  a  systematic  error  or  some  environmental  effect.  It 
is  normal  proceedure  before  every  cruise  to  calibrate  all 
temperature  sensors  in  the  laboratory  to  insure  that  they  read  the 
same.  This  was  done  for  MAGAT,  including  the  IR  thermometer.  A 
water  bath  with  an  immersed  platinum  thermometer  was  used  for  the 
IR  calibration.  The  two  temperatures  did  not  differ  by  more  than 
0. 3°C  from  0°C  to  40°C,  and  the  differences  could  be  accounted  for 
by  difficulties  In  mixing  the  water  to  ensure  the  bulk  and  surface 
were  in  equilibrium.  We  are  confident  that  any  difference  in  bulk 
and  IR  temperatures  measured  at  sea  are  not  instrument  problems. 

In  Figure  7,  we  show  the  measured  temperature  difference, 
IR-bulk ,  as  a  function  of  time.  Gaps  In  the  data  appear  for 
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times  when  the  ship  was  at  the  dock  or  the  measurement  system  was 
turned  off.  The  IR  temperature  is  always  lower  than  the  bulk, 
which  is  consistent  with  observations  that  we  have  made  on  all 
previous  cruises.  The  temperature  difference  varies  from  -0.6°C 
to  -3°C. 

We  have  checked  several  parameters  to  attempt  to  find  a  corre¬ 
lation  with  the  difference:  wind  speed,  swell  and  wave  heights, 
ship  location,  ship  speed,  bulk  temperature,  and  insolation.  No 
correlation  has  been  found  except  for  an  indirect  correlation  with 
insolation.  The  correlation  is  indirect  because  we  did  not 
measure  insolation  and  only  infer  it  from  time  of  day,  with  which 
the  IR  temperature  shows  a  correlation.  However,  the  possible 
insolation  correlation  is  opposite  to  what  should  occur  since 
absorption  of  solar  radiation  would  raise  the  skin  temperature 
with  respect  to  the  bulk. 

In  Fig.  7  we  also  show  a  rough  schematic  representation  of  the 
cloud  cover.  The  temperature  difference  is  lowest  when  the  skies 
are  overcast.  This  affect  can  be  accounted  for  by  reflection  of 
the  cloud  radiance  off  the  sea  surface  for  which  a  correction  can 
be  applied.  The  effect  of  such  a  correction  would  be  to  increase 
the  temperature  difference,  removing  some  of  the  fluctuations  but 
not  the  difference. 

Results  from  Wesely^-^)  measuring  IR  and  bulk  temperatures 
on  a  calm  cooling  pond  show  temperature  differences  from  0.3  to 
1.5°C  with  water  temperatures  varying  from  0  to  40°C.  The  effect 
should  be  smaller  for  a  wavy  surface.  Apparently,  we  cannot 
explain  our  results  on  the  basis  of  the  heat  transfer  rate 
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through  the  thermal  skin. 

As  of  this  time,  we  do  not  know  why  the  IR  temperature  is  as 
much  as  3°G  lower  than  the  bulk.  In  order  to  utilize  the  bulk 
model,  the  bulk  water  temperature  must  be  used  until  an 
appropriate  correction  to  make  to  radiation  temperatures  can  be 
found. 
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VII  CONCLUSIONS 


The  NPS  bulk  aerodynamic  model  for  calculating  the  optical 
index  of  refraction  structure  function,  Cj^,  works  quite  well. 

It  can  be  expected  to  predict  Cfj^  to  within  50 %  for 
homogeneous,  open  ocean  conditions.  In  coastal  areas  where  strong 
local  gradients  exist,  errors  as  large  as  a  factor  of  10  could 
result  if  meteorological  data  is  obtained  near  the  gradient. 
Weather  fronts  can  also  be  expected  to  produce  errors  but  there 
the  gradients  will  be  much  weaker. 

The  appropriate  sea  surface  temperature  to  use  in  the  bulk 
model  is  the  bulk  water  temperature  averaged  over  the  first  few 
inches  below  the  surface.  Temperature  measured  by  an  IR 
thermometer  cannot  be  used  directly  in  the  bulk  model.  The  small 
scale  turbulence  in  the  atmospheric  surface  layer,  which  is  in 
direct  response  to  air-sea  transfers  of  heat  and  momentum,  appears 
to  be  In  thermal  equilibrium  with  the  bulk  water  rather  than  the 
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